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Abstract 

We present an analytic solution of a charged black hole with hyperscaling violating factor 
in an Einstein-Maxwell-Dilaton model where the scalar potential is key to the existence of a 
solution. This solution provides a candidate gravitational description of theories with hyper- 
scaling violation at both finite temperature and finite charge density. Using this background 
we explore certain features of these theories via AdS/CFT correspondence. Finally, we dis- 
cuss embeddings based on the well-known sphere reductions of ten and eleven-dimensional 
supergravity. 



1 Introduction 



Gauge/gravity duality [T] may be thought of as a practical tool to study strongly-coupled systems 
near critical points where the system exhibits a scaling symmetry. Generally at a critical point 
the system may be described by a conformal field theory (CFT). From the gauge/gravity point of 
view this means that the gravitational theory is defined on a metric which is asymptotically locally 
Anti-de Sitter (AdS). 

On the other hand, in many physical systems critical points are governed by dynamical scalings 
in which, even though the system exhibits a scaling symmetry, space and time scale differently 
under this symmetry. A prototype example of such critical points is a Lifshitz fixed point where 
the system is spatially isotropic and scale invariant, though there is an anisotropic scaling in the 
time direction characterized by a dynamical exponent, z. More precisely the system is invariant 
under the following scale symmetry 

r -»■ \ z r, Xi -> Xxi, (1.1) 

where t is time and Xi are spatial coordinates of a d dimensional space (in our case the d dimensional 
space is R d ). 

In the framework of gauge/gravity duality it is then important to find gravitational theories 
that provide a gravity description of Lifshitz fixed points. Indeed, such a gravity dual for the 
Lifshitz fixed point was introduced in [2] where the authors considered a gravitational theory 
which admits a solution exhibiting the above scaling symmetry (see also [3] for earlier work on a 
geometry with Lifshitz scaling.). The corresponding metric may be written as 

ds 2 = -r 2z dt 2 + r 2 J2 dx 2 + (1.2) 

i=i r 

In terms of bottom-up model building, a natural set-up to consider is Einstein-Maxwell-Dilaton 
(EMD) theory, a setting in which numerous studies of Lifshitz-like black brane geometries have 
already appeared [4 - 14 . As is well appreciated at this stage, with the inclusion of dilaton and 



Abelian gauge fields, it is possible to find even more sophisticated metrics generalizing Lifshitz 
which, on top of an anisotropic scaling, have also an overall hyperscaling factor. More precisely, 
one may have a geometry in the form of [6] 



ds 2 = r=^ \-r 2z dt 2 + r 2 J2 + ^ j > (1-3) 

where the constants z and 9 are dynamical and hyperscaling violation exponents, respectively. This 
is the most general geometry which is spatially homogeneous and covariant under the following 
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scale transformations 



t -> X z t, r ->■ A V, -» Ax i; <is d+2 -> A*cis d+2 . (1.4) 

Note that with a non-zero 9, the distance is not invariant under the scaling which in the context 
of AdS / CFT indicates violations of hyperscaling in the dual field theory. More precisely, while in 
(d + l)-dimensional theories without hyperscaling, the entropy scales as T d l z with temperature^], 
in the present case where the metric exhibits hyperscaling, it scales as T^ d ~ e '' z 
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An interesting feature of the above metric is that for the special value of the hyperscaling 



violation exponent 9 — d — 1, the holographic entanglement entropy 23,24 exhibits a logarithmic 



violation of area law [17||25], indicating that the background (1.3) could provide a gravitational 
dual for a theory with an 0(N 2 ) Fermi surface, where N is the number of degrees of freedom, or 
alternatively the number of colors for an SU(N) gauge theory. 

These observations indicate that backgrounds whose asymptotic behavior coincides with the 
above metric may be of interest to condensed matter physics, thus making it natural to further 
explore gauge/gravity duality for these backgrounds, some of the aspects of which have already been 



touched upon in 18 19 . As the eventual goal of this program is to make contact with condensed 
matter, a natural generalisation involves backgrounds with finite charge density. Recently, we 
have witnessed the first construction of finite temperature and finite charge density solutions in 



a top-down construction with probe D-branes 20 . Here, as a further step in this direction, we 



identify what we believe to be the first back-reacted hyperscaling violating backgrounds with non- 
zero temperature and charge. This extends hyperscaling geometries in EMD theories at finite 



temperature 17 18 21 to finite charge density. 

This article is organized as follows. In the next section we will find hyperscaling violating 
charged black brane solutions in an EMD theory with a non-trivial potential for the dilaton. We 
observe that the presence and the form of the potential is crucial to have the solution. We shall also 
show that at zero temperature where the solution is extremal, the near-horizon geometry develops 
an AdS2 factor indicating that theories with hyperscaling violation at finite charge density and 
zero temperature may be described by an IR fixed point. In section three we study holographic 
entanglement entropy for a hyperscaling violating black brane solution, where we explore how the 
non-zero charge affects the entanglement entropy. In section four we probe the background by 
a charged fermion where we investigate the possibility of having an O(N ) Fermi surface in the 
model. In section five we study optical conductivity for theories with hyperscaling violation. The 
last section is devoted to discussions, in particular, those related to candidate embeddings. 

Note added: While we were writing up we became aware of further Lifshitz solutions with 
hyperscaling violation |22|. 



^or example, see recent discussion in 15 . 
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2 Gravity solution 



In this section we introduce our Einstein-Maxwell-Dilaton model and illustrate hyperscaling vio- 
lating charged black branes for this theory. For the asymptotically Lifshitz geometry an instructive 
analytic solution for charged black branes has been found (l4j. To find such a solution one needs 
at least two U(l) gauge fields coupled to a scalar field. The first U (1) together with the scalar field 
is necessary to generate an anisotropic scaling, while the second one is required to have a charged 
solution. It is also possible to find asymptotically charged black holes as well. For this purpose 



one should add another U(l) gauge field to be responsible for the compacted horizon 14 . 

To find hyperscaling violating charged black branes we will closely follow the approach of 
[14]. We note, however, that in order to get a non-trivial hyperscaling factor, besides the action 



considered in 14 , the model considered here has to have a non-trivial potential for the scalar field 
as well. In what follows we only consider charged black branes with one U(l) gauge field. We 
expect that it is straightforward to generalize for black holes with more gauge fields. 



To proceed we consider a minimal model as follows 

r 2 

S = I d^xy^g R - + V(<f>) - Ij^e^Ff 

J 1=1 



(2.1) 



where we have not confined ourselves to any particular dimensionality. For the potential of the 
scalar field, motivated by the typical exponential potentials of string theory, we will consider the 
following potential] 

V = Voe^. (2.2) 

Here Ai, A2, 7 and Vq are free parameters of the model. We will see that with this simple potential, 
the model does, indeed, admit solutions with hyperscaling violation. 
The equations of motion of the above action read 

*V + ^f-9„ = + I £ <** fa* pv - 9 ^F 2 ) , 



i=i 



dV(<f>) 



1 n 

+ V ^ (V=9~e Xi *Ft U ) = 0. (2.3) 



4 

i=i 



Let us consider the following ansatz for the metric, scalar and gauge field 

ds 2 = r 2a (-r 2z f{r)dt 2 + + r 2 dx 2 ) , = <f>(r), F t rt ± 0, (2.4) 

and assume that the other components of gauge fields are set to zero. 

2 Note that we do not expect to find this potential in string theory truncations of AdS vacua, but those of domain 
walls, i.e. the near-horizon of various branes. 
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From the Maxwell equations of motion, using the above ansatz, one finds 



Fi rt = e -^ r ^-d)+z-d-i p ^ ^ 
Then by combining the tt and rr components of the Einstein equation one has 

R\-K = -\g" W) 2 - (2.6) 

On the other hand, by calculating the Ricci tensor for our ansatz we get 

R\ - Rl = -d{a + 1)(q + z - l)r- 2a f{r). (2.7) 



Plugging this expression back into the equation (2.6) one can determine scalar field as follows 



e 4> _ e (po r y/2d(a+l)(a+z-l) __ e <f>Q r P _ (2.8) 

Note that in order to have a well defined solution one has to assume (a + l)(a + z — 1) > 0. Indeed 
it can be seen that this assumption is a consequence of the null energy condition. More precisely 



consider a null vector as 14 £ M = (y/g^, v^**, 0), then 



T^C ~ K -H$ = d{a + l){a + z- l)r~ 2c 7(r) > 0, (2.9) 

which in turn translates into (a + l)(a + z — 1) > 0. 

To find the metric one utilizes the other components of the Einstein equations. Indeed for our 
ansatz one has 

Rl = -( a + i) r -*(d+2)-z-d+i ( r d( a +i)+z f ^y ^ ( 2 .10) 

where, as in standard practice, a prime represents derivative with respect to r. Therefore the xx 
component of the Einstein equations of motion reads 

a(d+2)+z+d-l ( \r p ~l<\> 1 JL \ 

(r^ +1)+ V(r))' = (q + 1) f ^_ - ^^e-^- 2d(a+1) ) • (2.11) 



Using the solution of (2.8) one arrives at 



( r d(a+i)+* /(r) )' = !_ ^oe r _ i ^ e _ Ai ^ v 2 r _ 2d(a+1) _ A ^ ^ (2 12) 



(a + 1) \ d 2d 
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which can be integrated to find the function / of the metric as follows 
/( r ) = - mr -**-*-* + 



-dOL-Z-d , V e I 



d(a + l)(7/3 + a(d + 2) + z + d) 



2 r -2a(d-l)-/3\i-2d) 

_ n 2 p^ Xi<t>0 - (2 1 D 

j-f Pi 2d(a+l)(a(2-d) + z-d- p\iY 1 } 

where m is a constant which can be related to the mass of the black brane. 



With this expression, we have found all functions involved in the ansatz (2.4). The resulting 
solution has seven free parameters: a, z, m, 0q, pi, pi and the radius of the metric which has already 
been set to unity. On the other hand the action has four free parameters Ai, A2, 7 and Vq. Therefore 
four parameters of the solution can be fixed by the four parameters of the model yielding a solution 
with three parameters which correspond to the mass, charge and the value of the scalar field, let's 
say, at horizon. 

Note that so far we have not used all equations of motion. Indeed, in order, to find the 
parameters of the solution in terms of those in the action, we now utilize the remaining equations of 
motion. To proceed, we note that in order to maintain the asymptotic behavior of the hyperscaling 
violating metric, as it is evident from the expression of /, one must fix 7 as follows 

2ct 

7 = -j. (2-14) 

From this expression it is clear that to find a charged black brane (hole) with hyperscaling violating 
factor, it is crucial to have a non-trivial potential. 

From the equation of motion of the scalar field, we get 

4/3 M Vsfi*** = <**F? U - tt^tTt! (2.15) 

i=i ^ 



d(a + l) (3 ) u ^ 1 \ d(a + l) 



This equation can be solved for parameters Ai,A2 and p\. Actually this equation may be solved 
in different, but rather equivalent, ways. In particular one may find2 



/2d(q + l) 



A 2a{d-\) + 2d / 2(q + z-l) 2 = 2V (z - 1)£V 

l ~ y/2d(a + l)(a + z - 1)' 2 V d(a + l) ' Pl " da + d + z-1 

With these expressions, pi remains as an undetermined free parameter which, indeed, can be 
identified with the charge of the solution. The last parameter, Vq, may also be fixed by setting the 



3 It is worth mentioning that instead of finding the parameters of the solution in terms of those in action, we 
have fixed the action in terms of parameters of the solution. 
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constant term in the expression of / to one. Doing so, one arrives at 

V = e^K- 1 ^) (da + z + d - l)(da + z + d). (2.17) 

With this equation we have fixed all four free parameters, though there are still unused equations 
of motion. It is, then, important to check whether the other equations hold without imposing any 
further constraints on the parameters of the solution. In particular one of the non-trivial equations 
need to be checked is the tt component of the Einstein equations of motion. Indeed it is easy to 
see that this equation is also satisfied without imposing any further constraints. 



Therefore, to summarize, we note that the action (2.1) admits hyperscaling violating charged 



black brane solutions which can be recast to the following forrrj^] 



ds 2 = (-r 2z f(r)dt 2 + ^^ + r 2 dx 2 

V r 2 f(r) 

9(l-d)/d+d j 

F Xrt = ^2(z -l)(z + d- 0) e V^-on^-e/^° r d+z-e-i_ 



2 rt 



Q y / 2(d -e)(z-6 + d- 2)e~V^r *° r -i'+<^o-^ i 



4>o r yj2{d-e)(z-i-e/d) ^ ^2.18) 



witrS 2 

/(0 = l-^ + ra( ^- 1) - (2-19) 

This is indeed a charged black brane solution whose radius of horizion, r#, is obtained by setting 
/ = which leads to the following algebraic equation for r H : 

r 2( d+z -e-i) _ mr d+z-e-2 + Q 2 = Q (2 20) 

The corresponding Hawking temperature is also found 



_ (d + z-ey H f _ (d + z -e-2)Q 2 2{e - d . z+l) 
4tt 1 d + z-e Th 



(2.21) 



Therefore the extremal limit is given by 



2{d+z _e-i) _ (d + z-9 -2) 2 

Tli ~ d + z-e Q ' [Z - ZZ) 



4 Note that in order to follow the standard notation in the literture we set a = —6/d, where 8 > is the 
hyperscaling violation exponent. Observe also that this solution is not valid for 9 = d where a = — 1. 

5 As a consistency check , we note that for d + z — 9 = 2, this solution reduces to an uncharged solution obtained 
in (18 . 
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in which case the function / reads 



2jd + z-6-t) (r H \ d +*-\ d + z-9 (rH y{d+ z -e-i) 
/_i d + z-6-2 \r) + d + z-e-2\r ) ' ^ j 

In this limit it is easy to see that the solution near the horizon develops an AdS2 x R d ~ l geometry. 
More precisely, in the extremal case one can consider the following change of coordinates 

6T 2 T 

r ~ TH= ^ + z-e){d + z-e-i)C t= e^' (2 ' 24) 



Then the near-horizon limit is defined by the limit e — > for which the solution (2.18) reads 



2-2- / —dr 2 + d^ 2 

^ = Th " \{d + z-9){d + z-6-\)<? + 



A Xr 



9(l-d)/d+d , 

y/2(z - l)(z + d-9) eWd-BHz-i-e/d)™ i 



(d + z-6)(d + z-6-l) r Q H d - x C 

Qy/2{d-6){z-6 + d-2) e~V^T ^ 1 
2r ~ (d + z-0){d + z-6-l) r %+d-o-3 f ' 

= e^rt {d - e){Mld \ (2.25) 

which is AdS2 x R d with two charges. 

Following the general idea of gauge/gravity duality we would like to consider the solution (2.18) 
as a gravitational background dual to a theory with hyperscaling violation at finite temperature 
and finite charge density. We note that the gravitational description of theories with hyperscaling 



violation at finite temperature has already been studied in 18 . Actually the aim of the following 
sections is to further explore the holographic description of the theories with hyperscaling violation 
when the effect of charge has also been taken into account. 

As an observation, we note that due to non-zero charge, the black brane could be extremal 
with near horizon of AdS2 as we just demonstrated. It means that at low temperature the theory 
flows to an IR fixed point, where the theory may be governed by a two-dimensional CFT. 



3 Entanglement entropy 

The aim of this section is to study the holographic entanglement entropy of a hyperscaling vi- 
olating charged black brane by making use of the AdS/CFT correspondence. To compute the 
entanglement entropy via AdS/CFT correspondence one needs to minimize a surface in the bulk 
gravity. More precisely, given a gravitational theory with the bulk Newton's constant GV, the 
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holographic entanglement entropy is given by f23l|24 



Area(7 A ) 

Sa = 4GV ' (31) 

where 7^ is the minimal surface in the bulk whose boundary coincides with the boundary of the 
entangling region. 

In this section we will only consider the extremal case where the temperature is zero and 
therefore we will be able to study the effects of non-zero charge on the entanglement entropy. In 
the non-extremal case where we have non-zero temperature, the results essentially reduce to those 
in 18 . This is due to the fact that, in this case, the charge effects are subleading to the thermal 
behavior of the entanglement entropy. 

To proceed, let is consider a long strip in the dual theory given by 

1 I 

2 < x 2 < 2' < < L for % = 2,--- ,d. (3.2) 
The surface is defined by a hypersurface whose boundary coincides with the above strip and 



has a profile in the bulk of the metric (2.25) given by X\ = x(r). The induced metric on this 
hypersurface is 



7 2 2 
ds- m = p 



1 + x' (p) 2 ) dp 2 + dx 2 



(3.3) 



J{p) 

where we have made a change of coordinate, r = and prime now represents the derivative with 
respect to p. Therefore the area of the hypersurface reads 

A = L d - 1 [ dp p^^f-i + x' 2 . (3.4) 



As is well known, the above expression for the area may be treated as an action for a one dimen- 
sional mechanical system where the momentum conjugate of the field x is conserved. So that 



x po 



V 7 / 1 + x' 2 V P 



(3.5) 



with po representing the turning point where x'(Po) 00. It is then easy to find the width of the 
strip as follows 

i d - e r l/2 {i,P,/Pn) 



Jo v 7 ! " ¥ (d - 9) 

where 

f(t = 1 _ 2 ( rf + z ~ g ~ 1 ) ( Po\ d+Z ~\d+z-e , d + z-6 f po\ 2[d+z ^ 1] 2 {d+z-e-i) 
' Ph d + z-e-2 \p H ) ? d + z-e-2\p H ) ? 

(3.7) 
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The area is then given by 



A = L^p ^ 1 / 1 ? d f ~ 1/2 &P°/P * l. (3.8) 



Now the aim is to perform the above integrals and eliminate po to find the entanglement entropy 
as a function of the width of the strip t. Of course in general it cannot be done analytically, though 
in special limits it it possible to do so. 

Let's assume that the width of strip is very small, which in turn results in the turning point 
being much smaller than the radius of the horizon, i.e. po <C pn- For this case when 9 ^ d— 1, at 
leading order, the finite part of the entanglement entropy is 



d-lpO-d+l 

while for 9 = d — 1 the entanglement entropy exhibits a logarithmic violation of area law 



S finitc « Ld ~^~ d+1 (i + a ^ dj ) Q 4^i*+*-e + ai (z, d, 9)QH^ d+z - d -^) , (3.9) 



S « -jpr In - + c (z, d)Q^—— + Cl (d, z)Q 2 ——. (3.10) 



This is in fact a signature that the dual theory has an 0(N 2 ) Fermi surface which can be probed 
by the enanglement entropy 17 25 . Here a,i(z,d,9) and Ci(z,d) are numerical constants. 

On the other hand if we want to consider a strip with large width, one should let the hyper- 
surface extend all the way to the deep IR region, i.e. po ~ pn- In this case, at leading order, one 
arrives at 



L d-i, 



d-6 



^finite » -T^T- Q 1 +^~ 1 - (3.11) 

It is worth mentioning that if we compute the Hawking entropy of the extremal charged black 
brane the resultant entropy is 

Sh-^Q^H (3.12) 

which has the same behavior as the above entanglement entropy. Note that here Vol^ is the volume 
of d dimensional space generated by x. 



4 Charged Fermion probe 

In the previous section we have studied holographic entanglement entropy for theories with hy- 
perscaling violation. If the theory has an 0(N 2 ) Fermi surface, the holographic entanglement 
entropy is able to probe it. Indeed, in this case the entanglement entropy exhibits a logarithmic 
violation of the area law. Actually for the model we are studying when 9 = d — 1 with or without 
charge, we have a logarithmic violation of area law, showing that it has an 0(N 2 ) Fermi surface 
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(see also 17,25 



On the other hand when 9 7^ d — 1, the entanglement entropy has a power law behavior 
which indicates that the dual theory does not have an 0(N 2 ) Fermi surface. Nevertheless one 
may wonder whether the model still has an O(N ) Fermi surface. Of course in this case the 
entanglement entropy cannot probe the Fermi surface. In fact, to see whether the system exhibits 
an O(N ) Fermi surface one could probe the system by a charged fermion. Then to see a hint of 
a Fermi surface, one should look for a sharp behavior in the fermionic retarded Green's function 



at finite momentum and small frequencies 31 37] (for a review see |38|). 



In this section we will probe the extremal hyperscaling violating black brane geometry by a 
charged fermiorj^j This can be used to read off the retarded Green's function of a fermionic operator 
in the dual theory, which is a theory with hyperscaling violation at finite density. Using the pole 
structure of the corresponding retarded Green's function, we gain an insight into whether, or not, a 
Fermi surface exists. Since we will have to deal with fermions, for simplicity, we will only consider 
the four dimensional (d = 2) case. It is, of course, straightforward to generalize our analysis here 
to arbitrary dimensions. 



To proceed, let us consider a four dimensional charged Dirac fermion on the background (2.18) 
in the extremal limit 



bulk 



d x^/—gi^ 



V a Da-D a Y a 



m 



(4.1) 



Here Y a D a = (e p ) a Y»[d a + \{u pa ) a Y^ - iqA 2a ], with Y^ = \[Y»,Y% 

It is important to note that in order to find the equations of motion for this action, one needs 
to use the variational principle, which for a space with a boundary, always comes with a proper 
boundary condition. Of course, the boundary term is not necessarily unique and, indeed, for the 



above fermion action it was shown in 39 that there are numerous ways to make use of different 



boundary terms to make the variational principle well-defined. Of course, a given boundary term 
may break the symmetries of the model. In what follows we will consider the standard boundary 
condition (in the notation of (39) ). The equation of motion is 



e,) a Yt 



da + \(u} pa ) a Y p ° -iqA 2a 



m H> = 0. 



(4.2) 



Here the non-zero components of the vielbeins and spin connections for the metric in equation 

r e/2~z 



(2.18) are 



V7 



(4.3) 



5 Fermions in asymptotically Lifshitz geometries have been sudied in 26 -30 . 
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and 



„iz-e\i 



Sat, 



2-9 



[bJir 



(4.4) 



2riz-l-9)""' v ~"' u 2 
Then through the choice of ^ = (—h)~ 1 ^ 4 e~ luJt+lk/x ^lj(r), the above equation of motion reduces tc^] 



m 

„0/2 



^(r) = 0, 



(4.5) 



where /i 2 = QyJ 2 ^-e ^ e ^ ^~ e) ^° . Note that we have normalized the gauge field in such a way 
that it vanishes at the horizon. 

To find the retarded Green's function, we will utilize a numerical method to solve the equation 
of motion. To proceed, it is useful to consider the following representation for the four dimensional 
gamma matrices 



-a 3 



r* 



-a 



ia 1 
ia 1 



r 1 



-a 2 
a 2 



-ia 2 
io 2 



(4.6) 



Due to rotational symmetry in the spatial directions we may set fc 2 = 0. Then using the notation 



1> 



(4.7) 



the equation of motion (|4.5|) reduces to the following decoupled equations 

1 



1 



^-^(u, + ^ 2 (l 



for a = 1, 2. It is easy to see that 



9 mo 



r 



*« = 0, 



(4.8) 



$ a ~ a a 



for r — >■ oo. 



(4.9) 



Note that the asymptotic behavior of different components of the fermions are independent of m. 



This is in contrast to the Lifshitz case where we have 28 



for r — > oo. 



(4.10) 



Actually, this means that in our case we can use both standard and alternative quantization to 
read the retarded Green's function. In other words, one may consider either of a a or b a as the 
source and the other as the response. 



7 Here we have set the radius of the horizon to one. 
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To find the retarded Green's function, following [34] , it is useful to set xi — V'l/V^ an d Xi — 
1P3/1P4 where t/Vs are defined via $1 = (ipi,ip2),$2 — ("03; V^)- These parameters satisfy the 
following equations 



r 6»/2 A1 \ r zJl/2 ' r / Ai r zjl/2 r ' 

r/" 2 ^ + ^ " (^75 " 7) Xi = ^5 + ^, (4.H) 

where 

fi = w + ff ^(l--L). (4.12) 

Using these equations the retarded Green's function is essentially given in terms of functions 
G\(k,u) and G 2 {k,u) wher^] 

G a (k,u) — lim x a , for a = 1,2. (4-13) 



with the ingoing boundary condition imposed at the horizon, which in our notation is 34 

Xa I horizon ^- (4-14) 

More concretely, the corresponding retarded Green's function may be given by 



G 2 (k,u) 



Therefore the spectral function reads 



A(k,co) = hm^Gxik,^) + G 2 (k, cj)^. (4.16) 

Having found expressions for the retarded Green's function, it is an easy task to further determine 
its behavior as a function of k and u. We have illustrated the behavior in figure 1 for m = 0, qfi 2 = 

Vs. 

We have plotted the spectral function for 9 = l,z = 2 and 9 = 0A,z = 1.4. To compare the 
results with those in the RN AdS case studied in [34] , we have also plotted this case as well. Note 
that in all cases we have z — 9 = 1. Observe that since the retarded Green's function is an even 
function of k we only considered k > 0. 

One observes that the effect of increasing 9 or z, is to broaden the sharp peak representing the 
0(N°) Fermi surface. Moreover, as we see, even though for 9 = 0.4, z = 1.4 the entanglement 



^Here we set k\ = k. 
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I III 



Figure 1: The behavior of spectral functions the for hyperscaling violating geometry with 9 = 
l,z = 2 (left), and 9 = 0.4,2; = 1.4 (right). To compare the results with that of RN AdS black 
brane we have also plotted this case which in our notation corresponds to 9 = 0, z = 1 (see [34]). 



entropy does not exhibit a violation of area law, at the probe limit there is almost a sharp peak in 
the fermionic retarded Green's function which might indicate that the system has a Fermi surface. 

For z = 2 and 9 = 1, where the entanglement entropy exhibits a violation of area law indicating 
a Fermi surface at 0(N 2 ), we still have a peak in the spectral function around u ~ 0, though it is 
not as sharp as that in the RN AdS black brane. 



5 Optical conductivity 

In this section we study the optical conductivity of a finite density system whose gravity dual is 



given by (2.18). The local U(l) gauge field in the bulk is dual to a conserved current of a global 
U(l) symmetry in the dual theory. In terms of the conserved current, the optical conductivity can 
be obtained from the Kubo formula as follows 

<r((*>) = ^(^(^)^(-^))retarded, (5.1) 

where J{u>) is the conserved current of the U(l) global symmetry evaluated at zero spatial mo- 
mentum. The right hand side can be calculated using the AdS/CFT correspondence. Since we 
are working in the classical gravity regime, the leading contribution of the classical fluctuations of 
the gauge field is of 0(N 2 ) which is due to the black brane background. In particular, it is then 
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interesting to determine the optical conductivity for the case of 9 — d— 1 where we have an 0(N 2 ) 
Fermi surface. 

To proceed, we consider small fluctuations of the gauge field in the Xi-direction, 5A Xl = a±. In 
general, turning on small fluctuations of the gauge field causes a back-reaction on the other fields 
of the model, so we have to solve the equations of motion for all the fluctuations. More precisely 
one may have 



A 2fJ , -> A 2u 8 u0 + a 



2 /iW/iO 



l-i ' 



MNi 



-»■ </> + 



(5.2) 



Since we are interested in calculating the two point function, it is enough to expand the action 
up to quadratic level in the fluctuations. To do so, it is convenient to first use the gauge freedom 
and set h r M = a r = 0. Moreover at the level we are interested in we may set 5<p = 0. With 
these assumptions it is straightforward to compute the quadratic terms of the action from which 
we can read off the equations of motion for the fluctuations. For our purpose, we will focus on 



the zero momentum case which means we may set a± = a(r)e 



-iuit 



In this case the gauge field 



fluctuations mix only with the h tx component of the metric fluctuations. Therefore we find two 
coupled differential equations. Nevertheless one may eliminate h tx from the equations leading to 
the following differential equation for the gauge fluctuationp 



d r [V=ge x ^g rr g™a'(r)} + 



Q 2 g rr gtt 



ggx 



uj 2 y/^gg a g xx e x '^ ) a(r) = 



(5.3) 



where Q = \ / —ge X2 ^g tt g rr F2 r t- For the solution (2.18) the above equation read 
d r [r 3 ^ +d - e fa'(r)} 
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z-3+9-d 



2Q 2 (d-6)(z -6 + d-2)r 



Then for the special case of 9 = d — 1 the equation reduces to 



u} 2 r z-8+d-5 



f 



a(r 



0. 



(5.4) 



d r [r^- 2 fa'(r)\- 



2Q 2 {z-l)r 



z-A 



uj 2 r z 4 



/ 



air 



0, with/ 



r 2z ■ 
' H 



r,iY+ 1 , Q 



(5-5) 



It is interesting to note that the dependence on 9 and d drop out from the equation and it simply 
depends on a given z. 

Now the aim is to solve the above equation to find a solution with an ingoing boundary condition 
at the horizon. Then the corresponding retarded Green's function can be read from the asymptotic 
behavior of the solution near the boundary. Indeed one finds 



a ~ A + 



B 



r 3(z-l) ■ 



(5.6) 



J For more details see for example 42 and also Appendix 6.B of 



3 Note that in what follows, for simplicity, we have set (f>o = 



43 
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Figure 2: The real and imaginary parts of the conductivity versus u for a fixed temperature. To 
find this plot we have set Q = a/3, and r H = 1.25. 



for t — y oo and at the horizon we get 



a~(r-r h ) C'+D^-Q^-D. (5.7) 

Therefore with a proper choice of the solution near the horizon (ingoing solution) the retarded 
Green's function is given by 

(Jl(w) Jl(-W))retarded ~ ~7- (5-8) 

In order to find the retarded Green's function we utilize numerical methods. To proceed, we 



note that the equation (5.5) has four free parameters, Q, u,z and r#. In fact setting z — 2, and 
using "NDSolve" in Mathematica the real and imaginary parts of the conductivity as functions of 
frequency at a fixed temperature may be found as depicted in figure 2. 

The real part of the conductivity has the typical shape and from the imaginary part of the 
conductivity we observe that there is a delta function behavior in the real part. It is worth noting 
that the delta function behavior is an artifact of simplifications in the gravity calculation which 
should be compared with a sample without impurities. Adding impurities would broaden the delta 
function into a Drude peak. From the gravity point of view, this can be done by imposing a non- 
trivial boundary condition on the source of an operator at the boundary so that the translational 



invariance is broken. In this case, the Drude peak will appear from gravity calculations as well 44 
On the other hand, for a fixed u one can find the behavior of the conductivity as a function 

of temperature. In particular, setting u = 0.001, the real part of the conductivity can be found 

numerically as shown in figure 3. 

Again using numerical methods the best fit one finds for the real part of the conductivity is 

Re(a) « -0.22 + 0A2r H . (5.9) 
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1.5 - 



1.0 - 

Re(tr) 



0.5 




1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

Figure 3: The real part of the conductivity versus r# for a fixed frequency. To find this plot we 
have set Q = a/3, and to = 0.001. The numerical solution is shown by dots, while the best fit is 
drawn by a red line. 



Of course one may use the expression of the Hawking tempereture to write the result in terms of 
T. Indeed, doing so, one finds Re(cr) ~ T 1 / 2 . 



6 Discussions 

In this paper we have considered an EMD theory with a non-trivial potential for the dilaton. We 
have found an analytic solution representing a hyperscaling violating charged black brane solution. 
To do so, we have observed that the minimal model should contain at least two gauge fields, one 
to produce an anisotropic scaling and the other to correspond to a charge. It is known that to 
generate the required Lifshitz anisotropy, the first gauge field always diverges at the boundary, 
while the second falls off as a normal gauge field. Although in this paper we have only considered 
black branes with only one charge, we stress that generalizations to black holes with more charges 
are expected to be straightforward fli). 



Following the general thrust of the gauge/gravity philosophy, the resulting background (2.18) 
may provide gravity duals for strongly coupled theories with hyperscaling violation at finite charge 
density and temperature. It is then natural to explore different aspects of the theory by making 
use of gauge/gravity duality. Note that in order to study the effect of the charge, we have only 
considered the extremal case. Of course, one could have repeated the analysis for non-zero tem- 



perature, though the results at leading order would be the same as those obtained in 18 . This 
is because the effects associated to the charge are subleading. This can be seen, for example, in 
the expression of the temperature when expressed in terms of radius of the horizon. Moreover, the 



entropy of the black branes (|2.18|) at leading order in charge is given by 

(6.1) 



S ~ T 2 



(d - 9)(d+ z - 9 - 2) ^ 2 2(9+i-d-*) 

: Q J- 



z(d + z 
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An immediate observation is that in order to have positive specific heat, at leading order we find 



the constraint d> 9 18 , which was our assumption throughout the present paper. 



We have used the extremal solution (2.25) to compute the entanglement entropy for a strip and 
we have been able to obtain the leading order correction to the entanglement entropy due to non- 
zero charge. We have observed that, even with non-zero charge, when 9 = d— 1, the entanglement 
entropy, at leading order, exhibits a logarithmic violation of area law indicating that the dual 



theory may have an 0(N 2 ) Fermi surface 11 This, of course, would not exclude the possibility of 
having a Fermi surface of the order of one for arbitrary values of 9. 

Indeed, in order to study an order one Fermi surface we have also to probe the extremal 
background with a charged fermion. In the process, we have shown the retarded Green's function 
of the dual fermionic operator may have a pole in terms of a momentum for u ~ indicating 
that there is the Fermi surface. We note, however, that as one increases 9 or z, the sharp peak 
becomes smooth and eventually, for large enough 9 and/or z, the peak disappears, suggesting that 
the system does not have a Fermi surface. 

As we have already mentioned, at zero temperature the geometry will develop an AdS% geometry 
with a non-zero gauge field at the near-horizon limit. This indicates that at low energy the theory 
is described by a two-dimensional CFT, which is dual to the AdS2 geometry. However, due to the 
non-zero gauge field, the theory becomes unstable for sufficiently large charges. This is due to the 
fact that the scaling dimensions of operators in the dual theory become imaginary. If the operator 
is bosonic we would have bose condensation, while for a fermionic operator it may back-react to 



modify the geometry. Basically, the situation is very similar to that in a RN black hole 38 



It is also important to note that due to the hyperscaling violation and the behavior of the 



dilaton field at large distances, the background (1.3) cannot provide a dual description of a theory 



in all range of energies from UV to IR 40 41] (see also [18]) . In fact, gravity on the background 
(2.18) may be considered as an effective theory which is valid over an intermediate energy scale. 

What about string theory embeddings? We respond by echoing the sentiment of others by 
saying that it would be interesting to see if one can embed our solution to string theory. Here, to 
highlight that we do not expect this to be easy, we add some comments about embeddings based 
on the well-known sphere-reductions of ten and eleven-dimensional supergravity. In particular, 
we will focus on the truncations to U(l) 2 , U(l) 3 and U{1) A gauged supergravities without axions 
which exist in in seven, five and four-dimensions respectively (45] . Starting with d — 5 in our 
notation, neglecting the obvious drawback of it having a relatively more complicated potential, 
the U(l) 2 seven-dimensional gauged supergravity cannot be truncated further to a single scalar 
while retaining two independent gauge fields, so we can quickly eliminate this case. 

The potentials for d = 2, 3 are both of the form 

V{<j>) = Vxe^ + V 2 e i24> , (6.2) 



11 Another interesting observation is the charge dependence of the corrections is non-analytic. 
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and further truncations to a single scalar and two gauge fields can be found. So now we can return 
to section 2 and generalize the potential by adding the required extra exponential term. However, 
for these embeddings, if we require that one of the pi of our ansatz is a free parameter, such that 
it corresponds to the charge of a black hole, we find that the only solution is a = — 1, which we 
have already excluded from our solutions. 

So, the simplest truncations based on sphere-reductions do not allow for embeddings. Given 
that we initially only considered an exponential potential, it seems promising that one could embed 
this using a domain wall type set-up, since such simple potentials typically arise there. However, 
we still require the exact field content of a single scalar and two gauge fields. We leave further 
exploration to future work. 
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